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Abstract 

The neutral leptonic sector of the Standard Model presumably consists of three 

neutrinos with non-zero Majorana masses with properties further determined 

— . by three mixing angles and three CP-violating phases. We derive the general 

D . renormalization group equations for these physical parameters and apply them 






^ 



X 



to study the impact of radiative effects on neutrino physics. In particular, we 
examine the existing solutions to the solar and atmospheric neutrino problems, 
derive conclusions on their theoretical naturalness, and show how some of the 
5^ I measured neutrino parameters could be determined by purely radiative effects. 

For example, the mass splitting and mixing angle suggested by solar neutrino 
data could be entirely explained as a radiative effect if the small angle MSW so- 
lution is realized. On the other hand, the mass splitting required by atmospheric 
neutrino data is probably determined by unknown physics at a high energy scale. 
We also discuss the effect of non-zero CP-violating phases on radiative correc- 
tions. 
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1 Introduction 

There is mounting experimental evidence |]l],H,^ that flavour is not conserved in the flux 
of atmospheric and solar neutrinos. The most plausible, simplest and best motivated 
interpretation for this phenomenon is that interaction eigenstate neutrinos mix in a 
non trivial way into neutrino mass-eigenstates with different non-zero masses, leading 
to flavour oscillations 0]. The theoretical scenario that results most economical in 
accounting for the observed neutrino anomalie^ assumes that the known neutrinos 
of the Standard Model [ug, i/^, z/,-) acquire Majorana masses through a dimension-5 
operator 0, generated at some high energy scale A (the model can also be made 
supersymmetric) . 

The 3x3 neutrino mass matrix A^^ is diagonalized according to 

U^MyU = diag(mi, 7712, '"^3), (1) 

and we can choose rrii > 0. The MNS 0] unitary matrix U relates flavour and mass 
eigenstate neutrinos according to Ua = Uaii'i- The 'CKM' matrix, V , is defined through 

U = diag(e'°% e'°^ e'°0 ■ V ■ diag(e-''^/2^ e-'^'/^ 1), (2) 

and we use the standard parametrization 

V = i?23(^i) ■ diag(e-^^/^ 1, e''/^) ■ R^.iO^) ■ diag(e^^/^ 1, e'^^/^) . ^^^^^^^^ (3) 

where Rij{9k) is a rotation in the i-j plane by the mixing angle Ok-, that can be taken 
< 6'fc < 7r/2 without loss of generality. Explicitly, 

/ C2C3 C2S3 526"*'' \ 

V = -C1S3 - siS2C3e*^ C1C3 - siS2S3e*'' S1C2 . (4) 

\ S1S3 - CiS2C3e'^ -S1C3 - CiS2S3e*'' C1C2 / 

where Si = sin 9i, Ci = cos 9i. 

For later use, it is convenient to define the intermediate matrix W as 

W = V- diag(e-''^/2^ e-''^'/^ 1) . (5) 

The phases ae,afj,,ar, in eq. (^ are unphysical and may be rotated away by a redef- 
inition of the flavour-eigenstate neutrinos h'e,i'^,i'T, so that U coincides with W. It 
-'^ Leaving aside the, as yet unconfirnied, LSND anomaly ||]. 
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is useful, however, to keep these phases for the discussions of the next sections. The 
phases 6, and 0' are physical and responsible for CP violation in the leptonic sector 
(although only 6 can have an effect in neutrino oscillation experiments). The physical 
phases and 0' (usually kept out of the diagonalization matrix) are extracted from 
the equation 

V^M„V = diag(mi, m2, ms), (6) 

where mi = rriie^^, m2 = 17126^'^ , rh^ = m^. 

Let us now summarize the experimental information on the neutrino sector as set- 
ting the following bounds: 

From SK + CHOOZ data [jl|,|| we learn that 62 is small, with 

sin2^2 <0.1 (0.2), (7) 

at 90% (99%) C.L., or sin^ 2^2 < 0.36(0.64). The smallness of 62 implies that the 
oscillations of atmospheric and solar neutrinos are dominantly two-flavour oscillations, 
described by a single mixing angle 9i (precise analysis of the data must also include 
the subdominant effects due to a non-zero 6^2 0). 

In this approximation, oscillations of atmospheric neutrinos are dominantly z/^ — > z/^ 



with IC 



5 X 10"^ eV^ < Aml^^ < 10"^ eV^ , 

sin^ 29i > 0.82 . (8) 

Several parameter choices are possible to interpret the oscillations of solar neutrinos 

The large angle MSW solution (LAMSW), with 

10"^ eV^ < Am^^i < 2 x 10"^ eV^ 

0.5<sin2 2e3 < 1. (9) 

The small angle MSW solution (SAMSW), with 

3 X 10"^ eV^ < Am^^i < 10"^ eV^ 

2 X 10"^ < sin^ 2^3 < 2 X 10"^ (10) 



And the vacuum oscillation solution (VO), with 

5 X 10"^^ eV^ < Aml^i < 1.1 x 10"^° eV^ 

sin2 203 > 0.67. (11) 

Other relevant experimental information concerns the non-observation of neutrino- 
less double /3-decay, which requires the ee element of the Ai^ matrix to be bounded 

as 
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Mee=\mi clcje'^ + ma c^s^e^"^' + mg sj e'^"^ | < 0.2 eV. (12) 

In addition, Tritium /5-decay experiments indicate rrii < 2.5 eV for any mass eigenstate 
with a significant z/e component [0. Finally, no experimental information is available 
yet on the CP-violating phases. 

To write the previous bounds, we have followed the standard convention that or- 
ders the three mass eigenvalues with rris as the most split eigenvalue, i.e. |Ar7i2il < 
lAmlal, l^'TT'iil (where Amfj = mf — rrij) and identifies |Am|;^| = Am^^;, |Am|;^| = 
Am'^^^. Note that observations require |A?7i2i| <^ lAmlgl ~ l^'^iil (see however |T3[|). 



We do not adopt any particular ordering for 7711,1712- Nevertheless, note that although 
the interchange of mi, 1712 leaves sin^ 2^i and sin^ 2^2 unaffected, it fiips the sign of 



cos 2^3, which is important for the MSW effect \T^]. 

The relative wealth of experimental data has motivated efforts on the theoretical 
side to understand the possible patterns of neutrino masses and to follow the hints 
that such patterns offer on the fundamental symmetries underlying them. A more 
modest approach, perhaps more powerful at this given time, is to study Standard Model 
(SM) [or Minimal Supersymmetric Standard Model (MSSM)] radiative corrections (of 
well known origin) to neutrino parameters. Although neutrinos are famous for having 
extremely weak couplings to stable matter, this is not so for all the virtual particles 
that may appear in loops. Moreover, as we have seen, some of the mass splittings 
suggested by the data are very small and could be easily upset by radiative corrections 
of modest size. It turns out that, in many cases of interest, radiative effects have a 
very significant impact on neutrino physics. 



Many papers have dealt recently with these issues p!5| , p!6|p!7| , p!8| , p!9| , pO| , pT| , p^ , p^j . The 
most often considered scenario is to assume that a (Majorana) effective mass operator 
for the three light neutrinos is generated by some unspecified mechanism at a high 



energy scale, A. If, apart from this mass operator, the effective theory below A is just 
the SM or the MSSM, the lowest dimension operator of this kind is 

5C = --Kij{H ■ Li){H ■ Lj) + h.c, (13) 

where H is the SM Higgs (or the MSSM Higgs with the appropriate hypercharge) , Lj 
are the lepton doublets, and Kij is a (symmetric) matricial coupling. The neutrino 
mass matrix is then Aiu = kv'^/2, with v = 246 GeV (with this definition k, and Aii, 
obey the same RGE). Below the scale A, the most important radiative corrections to 
Aiu, which are proportional to log(A/M^) (where the Z boson mass, Mz, sets the low- 
energy scale of interest), can be most conveniently computed using renormalization 
group (RG) techniques. In short, one has to run Ai^ from A down to Mz with the 
appropriate RGE to obtain the radiatively corrected neutrino mass matrix. These 
corrections have interesting implications for the neutrino masses, the mixing angles 
and the CP-violating phases. 

As we will see, and has been discussed extensively in the literature, even if one starts 
with degenerate neutrino masses at the scale A, the relevant radiative corrections (with 
a typical size controlled by the tau Yukawa coupling squared and a loop suppression 
factor) tend to induce mass splittings that can be too large for some of the proposed 
solutions to the observed neutrino problems. They are particularly dangerous for the 
VO solution to the solar neutrino problem, and just about right for the SAMSW 
solution, which is certainly suggestive. They tend to be too small to explain the 
atmospheric mass splitting, which presumably requires to be generated by effects of 
the physics beyond A. One interesting possibility that, beginning with degenerate 
neutrinos, can produce radiatively a mass pattern in agreement with experiment is a 



see-saw scenario, for which the RG analysis was performed in [p]5| , [Lq] . 

Degenerate (or nearly degenerate) neutrino mass eigenvalues, besides being a very 
symmetrical initial condition (as befits the physics at a more fundamental scale) offer 
two important bonuses: the mass splittings at low-energy would be entirely due to 
radiative corrections and the mixing angles can evolve quickly to some particular values. 
The reason for the latter is the following. In the presence of degenerate mass eigenvalues 
there is an ambiguity in the choice of the eigenvectors which translates to the definition 
of the mixing angles. When a perturbation is added that removes that degeneracy, the 
ambiguity is resolved in a well defined way (which depends on the perturbation) and 



a particular configuration of mixing angles is cfiosen. Tliis is exactly what happens 
when a (sufficiently precise) initial mass degeneracy of two neutrino masses is lifted 
by radiative corrections. In RG language, the mixing angles can be driven quickly to 
infrared quasi-fixed points. Thus, RG effects could in principle be capable of explaining 
some of the measured values of the mixing angles. 



In this paper we extend previous work of us ]T5||T6|,|13 and others |[18| , p!9|j20| , |21| , |22| , |23 
along these lines. Assuming three flavours of light Majorana neutrinos, we derive the 
general renormalization group equations (RGEs) for their three masses, three mixing 
angles and three CP- violating phases (information alternatively encoded in the RGE for 
the mass eigenvalues and the complex mixing matrix U). This form of writing the RGEs 
represents an advantageous alternative to using the RGE for Aii, and, after obtaining 
the radiatively corrected mass matrix, extracting from it the physical parameters in 
whatever parametrization it is chosen. These two alternative methods can be described 
in short as 'diagonalize and run' versus 'run and diagonalize'. One advantage of the 
method presented here is that it avoids the proliferation of unphysical parameters 
(notice that TWi, has 12 degrees of freedom, but only 9 are physical), which allows to 
keep track of the physics in a more efficient way. Another one is that it permits to write 
the RGEs in a quite general form, without any reference to a particular scenario, such 
as the SM, the MSSM, see-saw, etc. This allows to appreciate interesting features, e.g. 
the mentioned presence of stable (pseudo infrared fixed-point) directions for mixing 
angles and CP phases, which are not consequence of a particular scenario. Finally, 
this method allows to determine the physical conditions for the validity of certain 
approximations, such as the two-flavour approximation, in a reliable way. 

In section 2, these general RG formulas are extracted from the most general form 
of the RGE for the neutrino mass matrix Ai^- These can be particularized to any 
model of interest. We give as particular examples the SM, the MSSM and a see-saw 
scenario. Some conclusions that hold in general are also discussed. In section 3, using 
the equations derived in the previous section, we study the implications of radiative 
corrections for neutrino physics in the SM case (or the MSSM) neglecting the effect of 
CP violating phases, so that the mixing matrix U is real. Section 4 extends this analysis 
to the more general case of non-zero CP-violating phases. Section 5 is devoted to the 
conclusions and the Appendix contains some technical details regarding the stability 
conditions for the mixing matrix U in the general case, and how they are reached. 



2 RGEs for physical parameters 

The energy-scale evolution of the 3x3 neutrino mass matrix A^,^ is generically described 
by a RGE of the form (t = log /x) : 

^ = -{kuM, + M.P + P'^M,), (14) 

(in particular A^^ = Ai^, is preserved). For interesting frameworks, such as the SM, the 
MSSM or the see-saw mechanism (supersymmetric or not), nu and P are known explic- 
itly p^j25|JT5[| . In (|T^), the term nuM.y gives a family-universal scaling of M.^ which 



does not affect its texture, while the non family- universal part (the most interesting 
effect) corresponds to the terms that involve the matrix P. 

The goal of this section is to extract the RGEs for the physical neutrino parameters: 
the mass eigenvalues, the mixing angles and the CP phases. We find convenient to work 
out the RGEs for the completely general case first {i.e. without specifying the form 
of Ku and P), and later specialize them for particular cases of interest. This allows 
to appreciate interesting features, e.g. the presence of stable (pseudo infrared fixed- 
point) directions for angles and phases, which are not a consequence of the particular 
framework chosen. 

Using the parametrization and conventions of sect. 1, one gets from eqs. (0, |T^), 
after some amusing algebra, the RGEs for the mass eigenvalues and the MNS matrix 

drrii 



dt 



We have defined 



-2miPii - miRe{Ku), (15) 



§^uT. ae) 



p=^U\P + P^)U, (17) 



while T is an anti-hermitian (so that the unitarity of U is preserved by the RG running) 
3x3 matrix with 

'mf + ni^APij + 2mimjP*, + iQi 



Ti- = - 
'^ (m? — mj) 



1-1 



VijRe{Pij) + i[Vij]-^Im{P,,) + iQ,,, i ^ j, 



where 



and 



Q = -'U\P - P^)U, 



Vi, = 



With the above definitions we can write 



rrii + irij 



rrii — rrii 



U^PU = P + iQ, 



(19) 
(20) 

(21) 



with P and Q hermitian. Note that the RGE for U does not depend on the universal 
factor Ku, as expected. 

From eqs. (p!6|-pOD, we can derive the general RGBs for the mixing angles, the 
CP-violating phases and the unphysical phases. For this task, it is useful to define 



T,i = T,ie^(<^-'^')/l n, = n,e'^/\ T,, = T,,/^'/^ 



) -^31 — -^31 



32 — -'326 



:22i 



which allows to conveniently absorb^ the phases (j) and <^' everywhere. Then we find 
the following expressions for the RGEs for the mixing angles: 



d9i _ 1 

dt C2 



Re 



S3T3I — C^Tr 



32 



(23) 



dt 



-Re 



(C3T31 + S3T32) e" 



■iS 



(24) 



dd, 
dt 



Re 



C2 



C2T21 + (S3T31 - C3T32 j S2e~ 



the RGEs for the CP-violating phases: 



dS , 

— - = Im 
dt 



1 ^ f V21 

T21 + ' 

S3C3 



\C1C2S2 



e-^ - ''^^ 



21 



S1C2S3 



)f„-( 



'3V2, 



22 



+ 



V, 



22 



\S1C2C3 C1C2S2 



Tc 



31 



= Im 

2 dt 



-'21 

■S3 



^3*2 C3K2* 



C1C2 S1C2S3 



21 \ T^ I ( "^31 
-'32 + 



V* V* 



2 dt 



Im 



^T2i +(^+ ^'' 
C3 



\C1C2 S1C2 



T32 + 



C1C2 S1C2 



21 I rp 
-'31 



+ Q33 - Q 



11) 



\C1C2 S1C2C3 



T- 



31 



+ Q33 - Q 



22i 



(25) 



(26) 



(27) 



(28) 



^The formulas we have written do not depend on our convention rrii ^ 0, which is just a definition 
for <j) and 0', and we could change at will the sign of rrii. The definition of the quantities Ty makes 
them independent of such choices of sign conventions. 



and the RGEs for the unphysical phases: 



~dt 



Im 



-T21 + ' =^2 



CsKl 



C3S3 



VciC2 



T32 + 



/K 



31 



S 3^2*2 



(ia„ 1 

— Im 



S1C2S3/ \C1C2 S1C2C3 

V;,f3i + 1/3*2X321 + ^33 . 



^31 



+ Q 



33, 



-Im 



(29) 

(30) 
(31) 



dt S1C2 

dOr 1 

dt C1C2 
Some comments are in order at this point: 1) the RGE (ITO) for U is only satisfied 

if the unphysical phases are included [that is, in the general case W does not satisfy 

eq. ([16|) ]. The reason for this is that the phases ae,a^,ar, that translate between U 

and the physical mixing matrix W, depend on the details in A1,y. When these details 

change through RG running the phases to be rotated away also change [thus eqs. (^9|- 

0)]. However, the explicit RGEs for the physical parameters rrii and 5,0,0' do not 

depend on ae, O/^, Ot as it should be the case; 2) from (p6|-^) we also see that, if 

no phases are present in U originally, they are not generated by radiative corrections 

(unless P contains phases. If it does, one could generate radiatively small non-zero 

phases in U); 3) from the structure of Tij (or Vij) we see that large renormalization 

effects can be expected in two cases: a) if some couplings in P are large (not the case 

of the SM) or b) if there are (quasi)-degenerate mass eigenstates at tree level (causing 

|V,,| > 1). 

An aspect of central importance in the discussion of the RG effects is that when the 

neutrino spectrum has (quasi)-degenerate eigenvalues, one expects large (even infinite 

for exact degeneracy) contributions to dU /dt. The reason is the following. When two 

mass eigenvalues, say i and j, are equal, there is an ambiguity in the choice of the 

associated eigenvectors, and thus in the definition of U . In particular, the columns 

i,j could be rotated at will, i.e. the matrix URij, where R is an arbitrary rotation 

in the i-j plane, will also diagonalize the initial Ai^, matrix. When the perturbation 

due to RG running is added, the degeneracy rrii = uij will be normally lifted and a 

particular rotation of the columns i, j will be singled out: U' = URij, giving T-- ^ 



and removing the singularity in eq. (pig). The form of R is thus determined by the 
matrix P. In the Appendix we give the derivation of R for a general P. The matrix U' 
may be very different from the form of U originally chosen, thus the initial big jump 
in the evolution of U. It is easy to see, however, that the subsequent evolution of the 
matrix U will be smooth, even in the presence of a large | Vjj|. If the initial degeneracy 
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is not exact, but nearly so, then |Vij| ^ 1 and the initial U, whatever it is, will be 
rapidly driven by the RG-running close to the stable U' form. Thus, U' plays the role of 
an infrared pseudo-fixed point. Therefore, if some initial neutrino masses are (exactly 
or approximately) degenerate, we have the interesting possibility of predicting some 
low-energy mixing angles or CP phases just from radiative corrections. This possibility 
will be exploited in the next sections. 

It is also useful to indicate that, in cases for which one particular Tij gives the dom- 
inant contribution to the RGEs, the following quantities are approximately conserved 
(here k can take any value): 

\Uk^\'+\UkJ\^ , (32) 

Im{UlU,,) , (33) 

[V,,]-ii?e(4-) , (34) 

as can re readily proved by using (^) and keeping only terms ~ Tij (for the last 
conserved quantity we assume that dP/dt is never of the order of Tj^, which is quite 
reasonable to expect). 

We turn now to particularly interesting scenarios. In the SM the expressions for 
nu and P are given by [p4| , ^ 

1 



Ku 



^[3^2_2A-6/i?-2Tr(YtYe)], (35) 



167r2' 

where (72, A, ht,Ye are the SU{2) gauge coupling, the quartic Higgs coupling, the top- 
Yukawa coupling and the matrix of Yukawa couplings for the charged leptons, respec- 
tively; and 

P = ^YtYe ^ ^diag(0, 0, 1), (36) 

where hr is the tau- Yukawa coupling. In the following, we will work in the approxima- 
tion of neglecting the electron and muon Yukawa couplings. 
In the MSSM one has instead 



/€[/ 



167r2 



U + Q9I - 6-%] , (37) 

5 sm p 

where gi is the t/(l) gauge coupling, tan/? is the supersymmetric parameter given by 
the ratio of the two Higgs vevs; and 

1 Y'I'Y 1 /?2 

P = -TTT^^^ ^ - 1^ 2 ;^ diag(0, 0, 1). (38) 

lovr^ cos^ p Idtt"' cos^ p 
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Note that we could have written these equations in terms of the supersymmetric cou- 
phngs ht = ht/sin(3, hr = hr/ cos P and Yg = Ye/ cos /9. 

In both cases, the previous formulas (|TB|, |TB[) apply with [see (P^-pB])] 



P — —K \TT |2 



-^ a — ^^: i 



T- 

-J- 1 



V 



rL-T- 



v.,Re{u;,u,,) + zV-'im{u;,u^ 



'3j, 



where 



327r2' 



for the SM case, and 



1 



hi 



rxjf 



(39) 
(40) 
(41) 

(42) 

ifi 2 2^' (43) 

lOTT^ COS"^P 

for the MSSM (note how t^j. is enhanced for large tan/5). For future use we also define 
the related quantity 

er = K Aog{k I Mz). (44) 

Hence, the RGEs for the neutrino masses given by eq. (pISj), as well as the RGEs for the 
mixing angles, the CP- violating phases and the unphysical phases given by eqs. (p3|-pl|), 
hold with the substitutions of eqs. (|39| , |iT| ). 

One can apply the general results to study other cases of interest beyond the SM 
or the MSSM. One particularly relevant case is the see-saw scenario |^. The model 
includes three heavy right-handed neutrinos Ni with Majorana masses given by a 3 x 3 
matrix M. with overall scale M ^ Mz such that the see-saw mechanism is imple- 
mented. Above M, the effective neutrino mass matrix is given by 



Mu = m^M ^rriD, 



(45) 



where 771^5 is an ordinary Dirac mass matrix coming from the conventional Yukawa 
couplings between the left-handed neutrinos, UgjU^jUT-, and the right-handed ones. The 
running of Ai and m£, is of the general form 

dM 



and 



dt 

drriD 
dt 



MPm + PIjM, 



moin'^h + Pd] 



(46) 



(47) 
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where k'u gives a family universal contribution. Combining (^) and ( ^7]) with ( ^Sf ) one 
gets an RGE for TWi, above the Majorana scale M of the form (0) with 

P = mB'PMm^-PD, (48) 

where we have assumed that mj:) has an inverse. Applying the general formulas to this 
particular case one can then easily reproduce previous results derived in the literature 
for this kind of scenarios IIT5|,|T6i . 



3 Implications (real case) 

The case in which CP-violating phases are zero, and thus the 'CKM' matrix V is real, 
is the one most extensively studied in the literature, perhaps because there is not any 
information yet about leptonic CP violation. This section is devoted to the study of 
this particularly interesting scenario, assuming that below the high-energy scale. A, the 
effective theory is just the SM or the MSSM, with an effective Majorana mass matrix 
for the neutrinos. It is convenient in this real case to work with the mass eigenvalues 
fhi defined in eq. (@) 

fhi = mie*'^, 1712 = m2G ? ^3 = i^s, (49) 

where 0, 0' = or vr. All m^'s are now real, but rhi and rh2 can be negative. We also 

define the quantities Vjj as 

~ rhj + fhi , , 

v., = ^ ^. (50) 

VTii — rrij 

Then, neglecting all the charged-lepton Yukawa couplings except hr, the RG equations 
for the mass eigenvalues and the matrix V read 

drhj 



dt 



2vmjV3j - rriiKu, (51) 



f^vr, ,5.) 



where Ku is given in eqs. {^^ ^7p, Kj- is given in eqs. ( ^^Hj ) and 



T^i = 0, 



Tij = KrVijVsiVsj (ij^j). (53) 
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The RGEs for the mixing angles are 

d9i 



, -/trCi(-S3V3iV31 +C3V32V32), (54) 

do 

—i = -K^CiC2(c3V'3iV3i + S3V32V32), (55) 

ijO 

-A = -K^{ciS2S3VsiV3i - C1S2C3V32V32 + V^3lV32V2l) • (56) 

dt 



From eqs. ( p^ p^D , the first conclusion is that the radiative corrections to the matrix 
V will be very small unless Ik^VjjI <; 1 for some i,j. This generically requires mass de- 
generacy (both in the absolute value and the sign), except for the supersymmetric case 
with very large tan/3, and thus large k^-. The interesting cases then are a completely 
degenerate spectrum ml ~ rrag — wig or an intermediate one with mf ^ m^ 9^ m'^. 

The initial degree of degeneracy of two given neutrino masses, say i,j, plays a 
crucial role for the subsequent RG evolution of V. In particular, it is important how 
small the initial mass splitting Amlj is compared with the typical one generated by 
RG evolution \6iiGE^i^ij\ ~ I^tI^^O) where mo is the average mass scale of the two 
neutrinos. There are four possibilities: 1) If Ami, — 0; there is exact degeneracy 
and V jumps immediately to V (note that we could have redefined the initial V to 
coincide exactly with the stable matrix V to start with, and then the evolution is 
smooth). 2) If lAr/ij-^- I < 2|kt-|?71o (i-e. |k^Vjj| ^ 1), V will quickly reach V (we will 
be more precise in the next subsection). 3) If 2|/€^|mo < lAm^yi < 2|eT-|mo {i.e. , 
[log(A/M^)]~^ < Ikt-VjjI ^ 1), V^ will tend to approach V but may not reach it if the 
total amount of running is not long enough (other possible effects in this case will be 
discussed later on). 4) Finally, if lAm^yi > 2\er\mo, no dramatic effects are expected 
from RG corrections to V (even if still the masses can be considered degenerate in the 
sense \Amlj I/ttiq <^ 1). 

For the discussion of cases of physical relevance, we will assume that the solar 
neutrino problem is solved by one of the standard solutions, so that |Am|^| -C jAmlgl ~ 
lArrig^l. In some of the cases that we will analyze (in fact, in the most interesting cases) 
the initial values for some of the m,i's are degenerate, so there exists an ambiguity in the 
labelling, which generically is removed after RG running. If the initial rrij's are only 
approximately degenerate (which we call quasi-degenerate), in principle there is no 
such ambiguity, but the RG running may alter the relative size of the Am? splittings, 
thus requiring a relabelling of the ms eigenvalue at low energy. Let us study in turn 
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the impact (and physical imphcations) of the RG running on the neutrino masses and 
the mixing angles. 

3.1 Mixing angles 

As we have discussed, in cases with some (sufficiently) degenerate neutrino masses, 
large RG corrections to the matrix U have the effect of driving it close to a stable form 
(an infrared pseudo-fixed point). This rises the interesting possibility of predicting (at 
least partially) low-energy mixing angles just from radiative effects. In the event of 
rrii ~ rrij, the stable form of V, say V, is characterized by the condition T-, ^ 0, which 
removes potential singularities in eq. (|5^, and requires VJj'j = or V^j = 0. The sign 
of the initial Kt-Vjj will determine which one is realized, as we will see shortly. 

We can estimate more precisely how V evolves to V in the following way. Suppose 
there are two nearly degenerate neutrinos with |eT-Vjj| ^ 1, so that this term dominates 
the RGEs of Vsj, V^j and Vjj itself: 



«j 



dt 



^r {V,] - Vi;) . (5^ 



Eq. ( p7D implies that, in the infrared, Vsj -^ (Vsj -^ 0) for positive (negative) initial 
KrVij. Without loss of generality, we can choose here the labels i,j so that vVj-j > 0, 
where the superscript (0) denotes the initial (high-energy) value. Hence, V^j -^ 0. 



Eqs. (|53, ^H) imply 

d 
di 



v., 



0. (59) 



The conserved quantity VsjVsj/Vij is the SM real version of the general conserved 
quantity given by eq. (Q). Notice that [Vj^]"^ = (m^ — m|)/(mj + mj)^ is a measure of 
the relative splitting of rh-j and rhj. For quasi-degenerate masses, [Vjj]^^ — Ainjj/Arn^. 
From (^) we find 

[v.,r^ = [vg^]-^-e.(^35-^3^), (60) 

where by {V^, — V^^ we mean the average value in the interval of running. When 
V ^>V' quickly, we can estimate that, with our choice of indexes 

K - ^l) - (^t- - ^1.) ~ -\ (61) 

13 



(note that V^j ^ by the running, while V'^^ ~ 1/2 is required to agree with experi- 
mental indications). 

Using (pup in the conservation equation (E^) we end up with the useful relation 






l-^.(^4-^3^)v!r ■ (62) 



This equation is interesting since it gives a measure of the change in VsjVsj, and thus 
in V, as a function of the initial Vjj and the amount of running. The stable form of 

V is achieved for VsiVsj -^ 0. From eqs. (|6T1 , |62|) we can evaluate the initial Vjj to get 

)/ 

'3j 



VsiVsj = Vgj VL /F, where F is an arbitrary factor 



^(o)t/(o)/i;. ^ ,^ f^(o) ^, 2(F-1) 



log M^ 

If V finishes close to the stable form, that means F ^ 1. For example, for F > 10 and 
A = 10^° GeV we need n^vf > 1. 



Notice that the condition to get F ^ 1, from eq. (|B5D, can be restated as [Am^J*^°) <^ 
2eT-?7iQ ~ 5RGE^f^%- This means that the final value for Am?- must be essentially given 
by 5RGE^rn'^ij, that is 

Am?- ~ 2KT-rn?\og——. (64) 

Mz 

This is interesting since it implies that if some |Vjj| is initially large enough to drive 
the matrix V into a stable form, not only the final mixing angles, but also the final 
A?7i?- splitting will be determined just by the radiative corrections. This result will be 
useful for the analysis of cases of physical interest. 

If V does not finish close to the stable form, iVsjVsjl may increase instead of going 
to zero (thus the value of F could be less than 1). For this to happen, one must start 
with HrV^i < KrV^j. Then [Vij]~^ decreases (as the mass eigenvalues get closer) and 
by eq. (|55|), iVsjVsjl initially grows. If the running stops before the eventual decreasing 



of iVsjVsjl below its initial value, then ^riV^j — V^^) > in eq. (|5^). Notice that a 
large increasing in VsjVsj needs a cancellation in that equation, namely ^I^t-I log 17- 



7(0) 



-1 



1-*^ 



(1 — F) Vij , which implies a certain amount of fine-tuning (the initial mass splitting 
must be of the order of the one generated by running, which implies a correlation 
between two quantities of totally different physical origins). 
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In the rest of the section we will consider the physics of the different types of 
spectrum which are relevant for the evolution of V (and thus of the mixing angles). 
They are the following: 

(i) ma^fhb^fhc^ \Vab\, iVfecl, iV^d > 1, 

(ii) -fha ~ -r/ifo cr^rhc^ \V ab\ > 1, 

{ill) rha^rhh^ ±7X1^ => \Vab\ > 1, 

[iv) —fha — rhh 9^ ±mc. (65) 

We leave the indices unspecified as they will be determined only after RG evolution 
(recall that m^ is defined as the most split eigenvalue). From the point of view of the 
running of the matrix V, what matters are the sizable |Vjj|. So, case (iv) is trivial (the 
matrix V changes little) and cases (ii)-(iii) can be analyzed simultaneously. 

Before embarking in the discussion of these cases, note that experimental observa- 
tions require 

V^i ~ (ss/v^, -ca/v^, ±1/72), (66) 

with different values for ^3 depending on the choice of solution to the solar neutrino 
problem. 

Case (i) 

In case (i), if |vVjj| > 1 for all z,j pairs, then V -^ V with V^jY-i^^ — > 0, and this is 
obviously incompatible with the desired form (|66D. Possible ways out would require 
starting with nearly degenerate masses but not so close as to drive V to its infrared 
fixed-point form. Such scenarios really belong in one of the next cases. 

Cases (a)- (Hi) 

Since m^ is by definition the most split eigenvalue, it is clear that in case (Hi) the label 
c corresponds to 3. This may not be so for case (ii), although we will see that a correct 
phenomenology requires it too. For the moment we will maintain the a, 6, c labels. 
In these cases the RGEs for the matrix V are dominated by the Tab terms. Explicitly, 

—— ^ 0. 67 

dt ^ ' 
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If initially KrVab > {KrVab < 0), then Vsb -^ (V^a -^ 0) in the infrared. From ( |66D 
we see that only V31 or V32 can be accepted to vanish [actually, this is guaranteed for 
case (in), since c = 3]. As we are free to interchange the labels 1,2, we choose V31 —>■ 1. 
On the other hand, the radiatively corrected masses are given by the expression (we 
absorb the universal scaling factor in a multiplicative redefinition of rhi) 

m,^m,{l + 2er{Vi;i), (68) 

where, as usual, {V^j) is the average value of V^^ in the interval of running. If |kt-Vjj| > 1 
then V -^ V quickly and we can set (^3^) ^ V'^^. 

Eq. (p6| ) and V31 — » imply that ¥% ~ V^'33 ~ 1/2. We have then, from eq. (p^). 



rrii -^ rrii, 

7712 -^ 777,2(1 + Cr), 

"^3 ^"^3(l + eT)- (69) 

The conventional labelling requires lAmlgl; l^^il ^ I ^^^21 1 5 so m^ and 777I must have 
an initial splitting large enough ~ Am'^^^ from the beginning (it will not be generated 
by the running)^. From eq. (|65|) we conclude that also in case (ii) c must be 3, in order 
to guarantee a correct phenomenology. On the other hand, as was discussed before, if 
initially |fi;rV2i| > 1, so that V reaches the stable form V, then the low-energy mass 
splitting is essentially the one generated by the running: 

A77721 ^ 2\er\ml. (70) 

Interestingly enough, this can be naturally of the right size for the SAMSW solution 
(3 X 10~^ eV^ < Am'j.^i < 10~^ eV^). For example, in the intermediate cases (in) this 
is achieved with sizeable values of the cut-off (A > 10^^ GeV) and/or working in the 
super symmetric scenario. 

The final (stable) form of V in this case is 



VI, 


VI, 


T/(0) 


vi. 


vi. 







vk 


T/(0) 

•^33 



V ^ Vi, ^'2 Vi^ , (71) 



3 



Conversely, if we do not start with such a splitting, eq. (69) indicates that mi becomes the most 



split eigenvalue, so it should be relabelled as 7713, leading to V33 -^ 0, which is not acceptable. 
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where V^j denote the initial Vij values. The values of the remaining VJ'- can be straight- 
forwardly written in terms of the former using the unitarity of the matrix V and the 
condition V^^ = 0. Moreover, we know that at low energy S2 ^ 0, si ~ ±l/v^, which 
allows to fill the third column of eq. (|7l|), and thus the complete final matrix V' 

/I \ 

\/' ~ 1/v^ 1/v^ , (72) 

VO Tl/v^ ±1/72/ 

This implies sin^ 203 ~ 0, which is also consistent with the SAMSW solution. It 
is worth- stressing that the value of sin^ 29^ is obtained just from the RG-running, 
independently of its initial valueQ. This was noted in ref. Jl^ for the (intermediate) 



case (in) , but clearly it also works for the completely degenerate case (ii). Let us also 
notice that the neutrinoless double /5-decay constraint implies that case (ii) is only 
consistent if m < 0.2 eV, see eq. (|12]). 

3.2 Neutrino masses 

The first consequence from eq. (|5TD is that neutrino mass differences get small modifi- 
cations, unless the scenario is supersymmetric and tan/? ^ 1, so that Hr ~ C^(l)- This 
means in particular that if the neutrino spectrum is hierarchical, i.e. ml < m\ <^ m|, 
satisfying Amfj ~ max{m^,m^}, then radiative effects are not going to appreciably 
change the spectrum and the mass splittings: \6rgeA^j\ *^ l^^^ljl- 

However, in a completely or partially degenerate neutrino scenario, i.e. ml ~ rn^ — 
m\ or m\ ~ m\ 9^ m| respectively, the shifts induced by the RGEs may be of the 
order or larger than the initial ones. They could also be larger than those required 
to explain solar neutrino oscillations. As for the mixing angles, the relevant cases in 
the study of the RGEs for the masses involve some complete or partial degeneracy. 
There is an important difference however. We learnt from the previous subsection 
that the mixing angles can only appreciably change if some Ikt-VjjI Z 1. For the mass 
splittings, there can be cases where, even for much smaller Ikt-Vj^I the RG effects are 
physically relevant. Actually, they can be very important too if the signs of the masses 
are opposite, rfij ~ —rhj, in clear contrast with the mixing angles (notice that in this 
case |Vjj| ^ 1). In particular, starting from a high-energy scale. A, down to Mz the 

^Note, however, that the value of V33 in this case is not affected by the RG running and thus has 
to be fixed by physics beyond A. 
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"solar" splitting Am2i can get a RG correction 

6RGEAml, ~ 4m^Kr{V^2 - V^) log -— , (73) 

Mz 

where (V3I — V^^) is to be understood as an average value in the interval of running. In 
degenerate scenarios this is typically much larger than the splitting required for the VO 
solution, unless (V^g^ — Vg^^^) ^ 0, i.e. (S1S3 — CiS2C3)^ ~ (—5103 — 0152^3)^ or, equivalently 



17,20 



cos^ 9i sin^ 62 — sin^ 9i ,_ ,. 

tan 29, ^ . , .%, '-■ 74 

sm 62 sm 2Ui 

This can be most naturally achieved with S2 — 0, sin^ 2^3 ~ 1, which amounts to a 
scenario close to bimaximal mixing, although it is not the only possibility. On the other 
hand, this way-out for the VO solution requires that V^2 ~ ^ii must be stable along 
the RG running. This implies in turn that both V32 and V31 must be stable. To see 
this notice that if V21 is dominant, then V§2 + V^^ ~ const., while if V32 is dominant, 
then V^i = 1 — V3I — V3I ~ const. As discussed in the previous subsection, and it is 
apparent from eq. (|52|), V will change little, unless some Ikt-Vj^I ^ 1, which corresponds 
precisely to a (at least partially) degenerate case in which the rhi and rfij masses have 
equal signs (|kt-Vjj| could also be sizeable if the scenario is supersymmetric with very 
large tan/5). Let us analyze all the possibilities. 

If fhi ^ 777.2, i-^- both have the same sign, then IV21I is the largest |V.jj|. IV21I 
should not be large enough to lead the matrix V to its stable form, as that would 
imply either V31 -^ or V32 -^ 0. But this means that Am'^i cannot be too small 
initially (or along the running). More precisely, it must be Arri^i ^ 2\er\'m?- In the 
case 777^ ~ 7772 ^ "^3 this is too large for the VO solution since irt^ must be at least 
~ A777^j^. The opposite case, with 777^ ~ m\ <C 7773, can be viable but only if the 
degenerate pair has very small masses, 777^ < 10~®eV^. 

With that caveat in mind, if 777^ ~ m\., the VO scenario requires 7771 ^ —7772, to 
avoid a disastrously large |V2i|. In that case, IV21I ^ 1 and V32 (or equivalently V31 
if it is 777i the one with the same sign as 7773) will be dominant. Still we need from 
eq. ([73|) (V^| — V^j^i) ~ with high accuracy, so one has to guarantee that the matrix 
V does not change much, to avoid dangerous changes in V^| — V^i- In particular, to 
avoid running into the stable form of V we require again A777I2 ~ 2m? Kr log j^. In the 
present case, this condition is easily satisfied for realistic A777I2, so it is natural to get a 



quite stable matrix V, as required. But still, V3I — V^i will (slightly) change according 
to 

^^i^i^ ^ 2K^Vi,Vi,Vs2. (75) 

Then, we can integrate eq. (|75D in leading log approximation and substitute in eq. ( [73D 
to obtain 

SrgeAttiI, ~ -4m\lVi^Vi2Vs2 log' — ~ --m\lVs2 log' — , (76) 

(we have used (^31^3!) ~ 1/16) which implies 

2—' 

atm 

This condition can be satisfied, but is barely compatible with a relevant cosmological 
role for neutrino masses and the VO solution (for the MSSM this conclusion is stronger 
as Kr is larger). On the other hand, condition (|77D is easily fulfilled in a partially 
degenerate scenario (the so-called pseudo-Dirac or intermediate case, m\ ^^ m^^ mg), 
where m? ~ Am^^^. Let us also remark that the splitting induced by the RGEs 
is potentially so large as compared to that required for the VO solution, that the 
constraint (^4]) must be satisfied with enormous accuracy [if some symmetry is invoked 
20| to justify ([74D, it must be either exact or minutely broken]. 

The previous paragraphs summarize and extend to the general case the results 



Am^, ^ AmL > 2-— ^«:^log' --. (77) 



from refs. |T5|JT6| , p!7| concerning the mass splittings induced by radiative corrections, 
especially in relation to the viability of the VO solution. 

3.3 The two-flavour case 



It has been claimed p!8|JT9| , |25[] that, working in the two fiavour (z/^, z/^) approximation, 
the RG running could drive the atmospheric angle sin^ 29i from a nearly minimal value 
at high energy to nearly maximal at low energy. This is an interesting possibility, since 
sin^ 29i is known to be nearly maximal. This subsection is devoted to this particular 
issue on the light of the previous discussions in this section. 

The analysis of refs. PB| , P^ , P5| was based on the RG equation for sin^ 29i [^| 



— sin' 29i = 2Kr sin' 29i cos' 29i , /^ ^ , '' . (7i 

dt ^ " ^ M33-M22 ^ 



19 



The observation was that if the diagonal elements of the mass matrix are close enough, 
sin^2^i could undergo a large increment. 

Now, starting from the general three flavour case, the two flavour approximation 
concerning RG running (not to be confused with the two-flavour approximation for 
oscillations) will be exact if 82 = 83 = 0, which is a stable condition. The matrix V is 
then simply given by 

V-\ ci si , (79) 

V -81 Ci / 

There is a physically interesting instance, namely the SAMSW solution, which essen- 
tially corresponds to this scenario. Then, from (p^, the RGE for 61 is 



—— = Kj. sm fc'i cos 61 V 32 = Kt- sm 61 cos 61 — ^ . ( 8U ) 

at m^ — 1712 

The corresponding RG equation for sin^ 29i is then given by 

— sin^ 201 = 2Kr sin^ 201 cos 261— . 81 

at rris — m2 

It is easy to check that eq. (|8l| ) is indeed equivalent to eq. ([78|) , but it is more practical. 



as it is written in terms of physical parameters. In any case, eq. (|80D is the more useful 
equation for the analysis^. Clearly, we can only expect a large modiflcation in the 
angle if |/€rV32| ^ 1. This requires 777,3 ~ ^2? unless the scenario is supersymmetric 
with very large tan/5 (^ 100). So, we must be in case (ii) (note that case (in) has c = 3 
by definition). However, our conclusions were that case (ii) can only be acceptable if 
labels a, b correspond to 1, 2, and thus c to 3, which is inconsistent with 7773 ~ 7772. The 
solution to this apparent contradiction is that the viability of this 7773 ~ 7772 scenario 
requires the running to stop before reaching the stability form of V, while in the analysis 
of subsection 3.1 we assumed that the stable form was reached. Let us see this in closer 
detail. 

If we start for example with ^i ~ (but not exactly 61 = 0, which is a stable value) 
and negative fi;^V32, the RGE ( pO]) will drive 61 — > 7i/2. Schematically, 



10^ j=, 72 - 1 



(821 



^ Incidentally, sin^ 26*1 — 1, i.e. maximal, is not a stable RG point, as has been argued in the 
literature |23] . From the general equations (pjjsfl) we can see that this is also the generic case in the 
3-flavour scenario, although one can easily determine particular conditions for stability. E.g. if V32 
is the dominant term and C3 = 0, then 9i is stable for any value. 
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(changing the sign of KTV32 would reverse the direction of the arrows). Notice that 
this is in agreement with our conclusions for case (ii): if |kt-V32| ~ 1, V33 -^ 0, which 
is inconsistent with experiment. However, there is an intermediate scale at which the 
mixing must be maximal, sin^ 29i = 1 (of course this does not occur if we start already 
near the stable 61 = 7r/2 value). This critical scale may be at low energy, and this is 
the possibility exploited in ref . |[l8|Jl9|j25[| . This of course requires some tuning of the 



initial parameters, and we will be more explicit about this shortly. 



Comparison of eq. (0) with the general equation for 61, eq. (54), shows that there 
could be other ways in which the two-flavour approximation for the running is correct. 
Namely, eq. (p^ reduces to ( pOD if one of the following possibilities occur 

(a) S3 = 0, 

(6) V31 ^ V32. 

These cases were not apparent at all from eq. ( |75D used in previous analyses. Now, 
scenario (a) is only acceptable if the S3 = condition is stable along the running, while 
Si is substantially changing. It can be checked from (|5^ -p6l) that this occurs only if 
S2 = 0. Therefore, scenario (a) can only work in the case previously commented, i.e. 
eq. (frOl). On the other hand, scenario (b) can work in many cases. 

Let us now estimate the amount of flne-tuning that these two-flavour scenarios 
require for the the atmospheric angle to be driven to maximal values at Mz thanks to 
the RG running. In particular, for case (a), we notice that the product V32V33 oc sin2^i 
must grow (V32V33 -^ V32V33/F) in a factor l/F, with F <^ 1. From the discussion in 
subsection 3.1, we know that this requires a suitable cancellation in eq. (|6^) between 
the initial mass splitting Am^ and that induced by the running 2eT-m^. Assuming for 
example a supersymmetric scenario (thus k^- < 0) with Vjj > 0, we flnd 

Am^ — 2e^m^ ^ , , 

< F < 1, (83) 



(we have made the estimate (V3I — V3I) ~ 1/2, which works well for small values of 
the initial mixing, as we have checked numerically). This shows that the flne-tuning 
between the two terms in the numerator of (^3]), which have completely different origins, 
is of one part in l/F. Therefore, the more important the increase in sin^ 29i, the greater 
the flne-tuning. Figure 1 illustrates well the behaviour of the running sin^ 29i and the 
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Figure 1: Running of sin 2^i in the two-flavour approximation. In the upper plot, we fix the initial 
value of the mixing angle but vary the initial mass splitting Am = 0.00393 eV in ±0.0005 eV. In the 
lower plot, we vary instead the initial value of the angle. We have chosen to = 1 eV and Kr = 6.2 x 10~^ 
to get maximal mixing near Mz- 
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fine-tuning involved. We show sin^ 26*1 as a function of tlie energy scale /i, from A = 10^^ 
GeV down to Mz (we let the angle run below Mz-, where it should really stop, just 
for illustrative purposes) according to eq. (PU]). The parameters have been chosen to 
reach maximal mixing near Mz starting from a small mixing angle at A. It is shown 
how, after reaching that maximal value, sin^ 26i goes down again, seeking the stable 
value [in accordance with the expected behaviour, see (|82|) ]. In the upper plot we 
demonstrate the dependence of the effect on the initial mass splitting: when it is varied 
slightly (leaving the initial angle fixed), the scale at which maximal mixing is achieved 
moves quite rapidly: an initial mass splitting slightly off and the mixing at Mz drops 
below the experimental lower bound. The lower plot shows instead the dependence 
with the initial condition for the mixing. When that initial value increases (decreases), 
the fine-tuning required is smaller, and this is reflected in a wider (thinner) half-width 
of the curve. 

Still within case (a), using again the conservation of V33V32/V32, we can easily 
evaluate the relation between the initial and final mass splittings as 

/S.m^{Mz) = Aml^^ ~ F Am'^{A). (84) 

Using eq. (^3|) we can obtain what is the value of Hr required as 

, , 1 ± F Am^ ^ ^ ^ 

For example, if F ~ 1/5, which implies an important increase in sin^ 29i with a 
moderate fine-tuning, and using 10~^ eV^ < m? < 4 eV^, Am^^^ = 10^'^ eV^ and 
10^ GeV < A < 10^^ GeV, we obtain 1.5 x 10"^ < {k^I < 0.28, which means that the 
scenario must be supersymmetric with tan/5 > 5. 

On the other hand, since V32 and V33 are not stable, we can evaluate from eq. ( [73l) 
the increment in the "solar" neutrino mass splitting Arri^i- Using now (V^l^^sl) ~ -^/2 
and (|85[) , we get 

SRGEAml, ~ AAml^ (86) 

which is too large. This is not surprising, as ml is stable but m"^ experiments a change 
of order Am^^^. Of course, one could still get Am2i of the right size tuning the value 
of nil so as to compensate the large RG correction (|86|) but that would be rather 
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unnatural. So, we conclude that the scenario (a), or equivalently eq. ([79|), cannot 
implement a substantial increase of sin^ 29i in practice. 

Let us consider now scenario (b). In principle, it can be acomplished in many cases. 
Namely, whenever ml < tti^ ^ fril (hierarchical spectrum); ml ~ TTig ^ rnl (inter- 
mediate or pseudo-Dirac spectrum); or ml ~ m,2 — mt\ (degenerate spectrum) with 
mi/m2 > and A?7i2i ^ ^^^32- ^^ the first two cases, V32 — —1. Therefore, if one 
demands a significant increase of sin^ 26^1 in eq. (|80D, then one needs |K^|log;|^ ^ 1, 
which signals the breakdown of perturbation theory. With regard to the degenerate 
case, since rhi and 777,2 have equal signs, the neutrinoless double /9-decay constraint im- 
plies that 777 < 0.2 eV, see eq. (|12]). This means that IV32I cannot be very large. Taking 
A777^^^ > 10~^ eV^, we get IV32I = 4?77^/A777^^^ < 160. Therefore, an appreciable vari- 
ation of sin^ 26'i requires Ivllog;^ ~ Am.'^^^/m? > 1/40, which means tan/? > 50, 
quite a large value. Anyway, in this degenerate case, IV21I ^ IV32I. Consequently, if 
IV32I is large enough to modify the V^i entries (which is our assumption), V21 will lead 
to V31 -^ (or equivalently V32 — ^ 0) rapidly, which requires an SAMSW scenario [see 
discussion after eq. ([71])]. But this is disastrous: from eq. ([7^), noting that (V3I — V^^) 
is necessarily sizeable, \6rge'A^i\ = C'(A777^^„), in disagreement with observations 
(again barring an unnatural cancellation between the initial and the RG-induced mass 
splitting). 

To summarize, the two-flavour approximation for the running of the 61 angle can 
be a good one if the previous (a) or (6) conditions are fulflUed. Then the size of Kr^32 
will determine whether sin^ 29i can get a substantial increase or not. If it can, then 
A77721 will get too large along the running. 

The possibility of obtaining an enhancement of sin^ 26i by the RG running seems 
so attractive that one may wonder if it could be realized in a 3-flavour scenario. Unfor- 
tunately, things do not seem to improve much. From the previous discussion, it is clear 
that the only possibility to consider here is 777^ ~ 7773 — 777^ (degenerate spectrum) 
with mi/m2 < and Ami^^ <^ A77732. Otherwise, the scenario is equivalent to the 
two-flavour case, which we know does not work, or it has too small Kt-Vjj to produce 
appreciable modiflcations in the angles. In the scenario selected only fi;T-V32 can be 
relevant. The effect of the running on the A77721 splitting is still given by eq. (0), 
which means that we will get |(5rg£;A7772i| = 0{Am,l^^) (unacceptable in any scenario) 
unless V^2~^ii ~ 0- This condition can be certainly arranged at low energy [i.e. at the 
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final point of the running) using S2 = 0, sin^ 263 = 1 (wliicli implies a nearly bimaximal 
scenario). However, if tlie running is going to modify appreciably the matrix V (which 
is mandatory to get an important enhancement of sin^26'i) this condition will not be 
stable. Notice that the dominance of Kt^32 imply that V31 ~ const., while Vs2,V3s 
should vary substantially. In consequence the final ^m\^ will be again too large. 

Schematically, the modification of the matrix U (starting with low sin^ 26i) will be 
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(This example has cos^i ~ 0, but the conclusions are the same for examples with 
sin Oi ~ 0). In the initial part of the running sin^ 26i experiments a large enhancement, 
so the parameters (kt-V32 and A) should be tuned to have the running finishing at the 
intermediate point. But in any case, it is apparent that V^l^^ai 7^ along the running, 
thus yielding too large Amii. 



4 Implications (general complex case) 

Non-zero CP-violating phases can have a non-trivial effect on the RG evolution of 
neutrino mixing angles and masses. Such effects have been considered previously in a 



2 generation case |21|. With the formalism presented in section 2 we can undertake 
easily the more realistic analysis for 3 generations. 

As in the previous section, we will assume that below the high-energy scale A, the 
effective theory is just the SM or the MSSM, plus the effective Majorana mass matrix 
for the neutrinos. As shown in sect. 2, the RGEs for the neutrino masses and the 



matrix U are given by eqs. (|T5|, ^6|), with the appropriate substitutions. Explicitly, 

drrij 



dt 



-2mi\Uzi\ 






rriiKu, 



^9) 



where ku gives a universal effect [see eqs. ( |55| , ^7\ )] and 



Tu = 



1 ij Kj 



V.,Re{U;,U,,) + iVTHm{U;,U,,)\ , {i ^ j). 



(90) 
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The explicit RGEs for the mixing angles, the CP-violating phases and the unphysical 
phases are given by eqs. (P5|-pi|), using everywhere the substitutions of eqs. ( P^ , |iT| ). 

As in the real case, large effects in the radiative corrections are expected when 
the neutrino spectrum has (quasi)-degenerate eigenvalues. When the degeneracy is 
sufficiently good (see previous section), U changes rapidly from its unperturbed t = 
form to a stable form that has non-singular dU/dt. To be more precise, if rrii ~ rrij, 
then |Vjj| ^ 1, and U is rapidly driven to a form that ensures Tij ^ 0. Since near 
t = 0, one has T^ ~ VijRe{U^^Uzj), the stable form, say f/', must satisfy 

R<U'lU'^,) = (91) 

for any pair z,j with mi ~ nij. Since the unphysical phase a,-, cancels out in the 
previous equation, we can also write 

Re{W'lW',,) = 0, (92) 

where the matrix W was defined in eq. (^. 

If all masses are quite different, then no dramatic RG effects appear and the discus- 
sion is similar to the corresponding one in the real case (see previous section). In the 
following subsections we then consider the two interesting cases of two-fold or three- 
fold degeneracy. The following general analysis contains as a particular case the real 
scenario considered in detail in the previous section. Note that now, two real scenar- 
ios with degenerate masses of either equal or opposite signs are treated together and 
correspond to choices of (p (or (p') equal to or tt, respectively. Note that, although 
the case with opposite signs has now a very large |Vjj|, having — 0' = vr gives a Tij 
under control, and the case is, of course, stable. 

As in the real case, for the discussion of cases of physical relevance, we will as- 
sume that the solar neutrino problem is solved by one of the standard solutions, with 

4.1 Three- fold degeneracy 

In the case mi ~ m2 — m^ ~ m, with |k^Vjj| ^ 1 for all 2,j, the matrix U (or 
equivalently W) quickly reaches a stable value U' {W') with 

Re{W'l^W'^2) = ReiW'l^W':,:,) = Re{W'l^W'^^) = 0. (93) 
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Using the fact that W'ss is real in the parametrization (||), the last two equalities in 
® imply 

iy^3 = or Re{W;,) = Re{W;,,) = 0. (94) 

If 14^33 7^ 0, eq. (|9^) and the first equality in (|93|) imply 



iy^i = 0, or iy;;2 = 0. (95) 

In other words, there must be some W^^ = (the ambiguity in the labelling of the 
three original eigenvalues is reflected in the ambiguity in the z-label of W^^ = 0). Then, 
by unitarity, the values of iM^^jP = |f^3jP are {0,x, 1 — x}, with < x < 1 (the value 
of X is determined by the initial form of the matrix W). So, from the expressions for 
drrii/dt we obtain the following low-energy shifts in the mass eigenvalues: 

Amj ~ me^{0, X, 1 — x}, (96) 

where m is essentially the initial value of the three masses at the scale A. 

Therefore, for a given value of x we know how the degeneracy in the mass spectrum 
is lifted and then the labelling of the three eigenvalues is determined according to our 
convention (thus U' is fixed). It is then easy to show that no value of x gives a stable 
U' in agreement with the available experimental information on the neutrino mixing 
angles. For x < 0.25 (or > 0.75) one has U^^ -^ and jf/sap > 0.75, in conflict with the 
limit cfc2 < 0.71 which follows from the experimental bound (^). For 0.25 < x < 0.75, 
f/33 -^ instead, and this is also incompatible with experimental limits: it implies 
either c^ — ^ [in conflict with ([^)], or c^ ^ [in conflict with (^)]. 

In summary, the case of three-fold degeneracy of initial neutrino masses leads to a 
stable U' which does not accomodate the values of 9i and Arnfj suggested by exper- 
iment. Note that this analysis contains and expands the subcases (i)-(ii) for real U 
analyzed in the previous section. This negative result has two obvious way-outs. If the 
spectrum is almost degenerate, but all |fi;,-Vjj| <^ 1, then the matrix U is very slightly 
renormalized, so it never reaches the stable (disastrous) form. In this case the RG cor- 
rections to mixing angles and phases are basically irrelevant. This extreme situation 
is difficult to implement in the VO scenario, since |fi;T-V2i| is necessarily sizeable. The 
second way-out is that |kt-V32| ^ 1 (even though IV32I ^ 1), while |fi;T-V2i| can be sig- 
nificant. This will occur if the Am32 splitting is arranged to sensible ("atmospheric") 
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values from the beginning. Then, the analysis for the mixing angles and CP phases is 
identical to a case with two-fold degeneracy, which is analyzed in the next subsection. 

4.2 Two-fold degeneracy 

We assume here that the initial Avrigg, A?7i|^ splittings are large enough to keep ms 
as the most split eigenvalue after the RG effects, while mi ~ m2- In this case, only 
|'«rV2i| ~ 1. This drives the original matrix U to a. stable form U' given by (see 
Appendix) 

U' = URu{T), (97) 

where i?i2(r) is a rotation in the 1-2 plane by an angle F, to be computed below. We 
see that the third column of U is not changed appreciably, in accordance with the 
fact that dUis/dt does not depend on V2i. The rotation angle F is determined by the 
condition (pi|) 

ReiU'*M,) = 0, (98) 

which gives 

This allows to determine U' completely and thus the physical parameters as a function 
of the initial conditions. Let us examine this issue in closer detail. 

With regard to the mixing angles, from the RGEs (^^- |25|) one immediately sees 
that 6i and 62 are not renormalized strongly {d6i^2/dt do not depend on T12) while 63 
is driven to a stable value, 6^ . More precisely. 



e'i^ ~ 62 

sin^ 29i^^ = p^ + q' , (100) 



where 

p = sin2Fsin[((/)-(/)')/2], (101) 

and 

q = sin 2F cos 2^3 cos[(0 - (f)')/2] + sin 2^3 cos2F. (102) 
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The CP- violating phases are also driven to particular values 5^-^^ , (p^^^ , 0'*^-^^ which are 
complicated functions of the initial values of the parameters in U: 

tan^(/) = ^^""^-r (103) 

q + p tan o 

tanfrc^(^) - 6'^f^)/2] - sm2e,sm[{(j) - (j)')/2] 

lan[[(p (p J/^J-^os2rsin2^3Cos[(0-0')/2] + sin2rcos2^3 ' ^ ^ 

and 

0(/) + ^'(f) _ 25(/) = (f) + (j)' - 2S. (105) 

The last equation follows from d{(f) + </)' — 26) /dt ^ 0, which can be checked directly 
from (|6|)-(^. 

As explained in section 2, see (^31), other quantities of interest which are con- 



served in this case are Im{Ul]Uk2)- For A; = 1, we obtain the conserved quantity 

sin2^3sin[(0-0')/2], (106) 



and, for /c = 2, 3, 



sm 



0')/2] cos (5 cos 2^3 - sin(5cos[(0 - 0')/2]- (107) 



These expressions can be sometimes conveniently used instead of the explicit formulas 

(iii-(ra- 



For physical applications it is interesting to use for 6i and 62 the values suggested by 
experiment]^, sin^ 26i ~ 1 and ^2 ~ 0, which permits to obtain more explicit expressions 
for 6*3 and the CP-violating phases at low energy. In doing this, it is important to keep 

corrections due to a small (but non-zero) value of ^2, as they can become dominant for 

if) 
some choices of the parameters. For the angle 6)^ , we obtain the remarkably simple 

expression 

sin^ 2^^^ = sin^ 2^3 sin^cf) - 00/2] + O(s^), (108) 

i f) 
which gives the final ^3 as a function of the initial 63 and (p ^ 4''- We see that, for 

the SAMSW solution of the solar neutrino problem, we must have either sin^ 2^3 ~ 

or sin^[(0 — 4>')/2] ~ at the scale A. On the other hand, if the solar mixing is also 

^Due to the stability of 0i and 62, this initial condition at A is approximately maintained by RG 
evolution down to Mz- The origin of such initial condition could be some flavour symmetry or perhaps 
a fixed point in the RG running above A. 
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maximal, sin^ 26^3 -^ ~ 1, then both sin^ 26*3 and sin^[(0 — 0')/2] must be of order 1 
originally. 

A word of caution should be said about the higher order term in (p,08|). The 



neglected term, calculated exactly, is [we introduce the short-hand notation ys = (0 — 
00/2 and a = sin 2^1 = ±1]: 

4^27^^^' (109) 



with 



A^i = c^cs + s^S5C2e3, 

-Di = -C<^C2S2e3 + 2(TS2(c<^C5C293 + S^Ss)-, 

D2 = C2C293 + 2crs2CsS2e3 , (110) 

where all the quantities are understood at the initial (high energy) scale. For generic 
values of the parameters, the term (|109| ) will be negligible as it is formally 0(^2). 



However, for some particular choices of the parameters the denominator in (|109|) can be 
also 0(^2) resulting in a non-suppressed correction in (|108|) . This occurs, for example, 
if 026)3 — ^^^ c^ ~ simultaneously. 

Expressions ( [108| ) and ( |109| ) for sin^ 2^3 , in conjunction with the conserved quan- 
tities ( |105| - |107| ), are enough to determine exactly the rest of parameters at low-energy. 
For example, we obtain 

tan-|(^'"-^-"')/2|^ "°^^^^^^<g + ^^' . (Ill) 

Note the l/s^ dependence that would generically drive (p^-^^ — (p'^-^^ -^ ±n. However, 
this is not necessarily the case because the smallness of S2 can be compensated by a 
small numerator in ( p.ll|) , as commented above. 

If = 0' initially, it is easy to see that the only physical parameter that undergoes a 
sizable change is ^3, with 9^ = ^3-l-r. This follows directly from U' = URi2{T) because 
now -Ri2(r) commutes with diag(e*'^, e*"^, 1). Moreover, eq. (|108| ) gives ^3 -^ C>{s2) ~ 
irrespective of the initial value of ^3. This is just right for the SAMSW solution of the 
solar neutrino problem and reproduces the result discussed in the previous section for 
the real case. We see however that this appealing possibility dissappears for generic 
values of and 0'. 

If — 0' = ivr instead, one gets tan2r ~ S2 sin 6. For ^2 = or 5 = this implies 
a stable matrix U from the beginning, see eq. (0). This is in correspondence with 
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the real case = 0, 0' = tt discussed in the previous section. For non-zero values of 5 
and/or ^2, however, the matrix U is no longer stable, although typically F, and thus 
the corrections to f/, will be small. We conclude that, the presence of a non-zero phase 
(5 = tends to destabilize this scenario, but the effect is small for small 62. 

Another interesting result that we can extract from eq. ( |103| ) is the possiblility of 
generating a non-zero phase 5'^^^ even if 5 = originally, provided the phases 0, 0' are 
non zero and satisfy sin[2(0 — 0')] 7^ 0. The resulting 5'^^\ when 5 = 0, is given by 



tan(5(^) = -^, (112) 



or, explicitly. 



^ r(f) _ -s^c^jcl sin 26^3 - 2(JS2 cos 29^^) 

~ Cis2 sin 2^3 cos 2^3 + 2(752(1 -s^cos^ 2^3) ^ ' 

This shows that one can generate radiatively even a maximal value for the CP- violating 
phase. 

4.3 Neutrino masses 

In the complex case, the running of the neutrino mass eigenvalues is governed by 

dm.; 



dt 



-2Krmi\U^i\^ - Kuiiii, (114) 



and many of the generic results obtained in the real case go over to the complex one. 
Some of the significant differences are discussed in this subsection. 

From (|114D we see that the phases and 0' do not affect directly the running of the 



masses. However, they have an important influence in the stability of a given choice of 
the matrix f/, and this in turn will influence the evolution of the masses. For example, 
in the degenerate case mi ~ m2 = m,Q, the initial choices — 0' = and — 0' = vr 
(which correspond to rhi ^ ±^2? already analyzed in the real case) give identical 
dnti/dt at t = 0. However, in the first case U changes abruptly, while in the second 
it evolves smoothly, and this has an evident effect on the subsequent running of mi 2. 
The phase 5, on the other hand, enters f/31 and f/32, and thus has a direct effect on the 
size of dm^i/dt, although its importance is somewhat reduced by the smallness of 62. 

One important instance in which the stability of U is crucial for the viability of 
the scenario is the VO solution to the solar neutrino problem. As discussed in the real 
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Figure 2: Allowed regions (in which Aml^ < 1.1 x 10"^" cV^) in the plane (sin^ 2^2, sin^ 2^3) for the 
SM case for sin^ 20i — 1 (upper plot) and 0.82 (lower) and a cut-off scale A = 10^^ GeV. The different 
strips correspond to different values of the phase S, from (darker color) to 7r/2 (lighter). 
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case, in a scenario with total or partial degeneracy of neutrino masses, the final Arrigi 
has to be much smaller than the typical size of the RG shift in the masses squared 
~ 2mQer- Hence, starting with rrii ^ 1712 ^ rris at the high energy scale A, to keep 
Amli under control requires an exquisite cancellation in the running of mi and m2. 
From (|114|) , this cancellation is 

mkf'-\uL\')\f.l^^. (115) 

In this condition, we have written U' and not U, because if the degeneracy is good 
enough U jumps to U' almost immediately away from A and remains close to it during 
the rest of the running down to Mz- By definition, U' satisfies the stability condition 

ReiU'l,U',2) = 0. (116) 

From ( |115| ) we find that, to keep the solar mass splitting under control, the mixing 
angles (at low energy) must be correlated according to 

cos^ 01 sin^ 02 - sin^ Oi 

tan203~ ■ n ■ oa ^- ^^^^) 

sm U2 sm 2&i cos 

This equation generalizes ( [7^ to the complex case. Eqs. ( |116| , |117D must be satisfied 
simultaneously for the viability of the VO scenario in totally or partially degenerate 
scenarios. 

From eq. (|117|) we can see that a non-zero 6 will infiuence the correlations that must 
occur. The effect is discussed in figure 1, where, for two different values of sin^ 29i (the 
maximal one, 1, and the experimental lower bound, 0.82) we plot the region of the 
plane (sin^ 2^2, sin^ 26^3) where the cancellation (|115| ) takes place. The different stripsf] 
correspond to different values of cos 5, from 1 (darker color) to (lighter), in 1/4 
steps. We set the high energy scale A at the typical see-saw value 10^^ GeV and require 
A?7i^2 ^ 1-1 ^ 10^"*^° eV^ with the initial neutrino mass itIq = Am'^^^ = 5 x 10~^eV^. 



For much smaller A the strips in figure 1 get thicker [it is easier to satisfy ( |115| ) because 



the size of radiative corrections is smaller] . Note that 6 = 11/2 is special because, in 
such case, eq. ( |117|) has in principle two solutions: a) sin^ 26'3 = 1 and b) s| = tan^ 9i. 
The region satisfying ( p,15| ) will consist of two strips centered around those two lines. 



^For a fixed value of sin 26*2 and sin 29i there are in principle 4 values of sin 263 satisfying (117), 
corresponding to the fact that the interchange cos 9 -^ sin 9 leaves sin^ 29 invariant. We do not plot 
the two solutions with S2 > C2, which would be in conflict with CHOOZ+SK data. Note that, for 
sin 201 = 1 the two remaining solutions merge in one. 
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In the case sin^ 26*1 = 1, the solution b) gives Sj = 1, which is beyond experimental 
bounds and only solution a), sin^ 2^3 = 1, remains (the region around that line has 
zero width and cannot be seen in the plot). For sin^ 26*3 < 1, solution b) can have 
S2 < C2 and we include it in the figure. The corresponding region is T-shaped like the 
one shown in figure 1, lower plot. 



On the other hand, the condition (|116|) for stability can be satisfied for any choice 



of mixing angles and 6 by adjusting — 0' to the apropriate value, which is given by 

, sin2^3(sin^ 6^1 — cos^6'i sin^ 6^2) — sin26'i cos26'3sin6'2 cos(5 

tan - /2 = ■ 00 ■ a ■ Ji ' ^^^ 

sm 2&i sm 62 sm 

This would be the complex version of the requirement of opposite signs for the two 



degenerate eigenvalues in the real case. Using ( |11?| ), the previous result simplifies 

further to 

tan[(0 - 0')/2] = --. (119) 

^^ ^' ^ tan (5 cos 26*3 ^ ^ 



Unlike the real case, to impose condition ( |118| ) for generic values of the parameters 



seems hard to justify from some underlying symmetry. 

5 Conclusions 

Assuming three fiavours of light Majorana neutrinos, with no reference to any particular 
scenario, we have derived the general renormalization group equations (RGBs) for 
the physical neutrino parameters: three masses, three mixing angles and three CP- 
violating phases [information alternatively encoded in the RGE for the masses and 
the complex mixing (MNS) matrix U]. This form of writing the RGEs represents 
an advantageous alternative to using the RGE for Ai,^. It avoids the proliferation of 
unphysical parameters, which allows to keep track of the physics in a more efficient 
way. It also permits to appreciate interesting features, e.g. presence of stable (pseudo 
infrared fixed-point) directions for mixing angles and phases, which are not consequence 
of a particular scenario. 

We have then particularized the RGEs for relevant scenarios. Namely, when the 
effective theory below A (the scale at which the neutrino mass operator is generated), 
is given by the SM or MSSM; or when Ai^, is generated by a see-saw mechanism. In the 
first two scenarios we have analyzed in detail the physical implications of the RGEs, 
separating the case where U is real {i.e. no CP phases) and the general complex case. 
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For the real case, we have noticed that if one starts with two masses suficiently 
degenerate (in absolute value and sign), say rrii ~ rrij, the U matrix is driven to 
a stable (infrared pseudo-fixed point) form, providing net predictions for the mixing 
angles. Whenever this happens the corresponding mass splitting. Am?-, is entirely 
determined at low energy by the RG running. This amounts to a very predictive 
scenario. Depending on what are the initial (quasi) degenerate neutrinos, the scenario 
can be realistic or not. In particular, starting with mi ~ m2 and an initial Am'^2 ~ 
Am^j^, we finish at low energy with a small "solar" 9^ angle and a Amgi sphtting 
just of the right size for the SAMSW solution to the solar neutrino problem, which 
is certainly remarkable. On the other hand, the radiative corrections to the mass 
splittings are potentially dangerous for the VO scenario, which becomes unviable if 
mf ~ 1712 ^ "^3 5 unless mi ~ —m2, the "solar" ^3 angle is close to maximal and the 
common mass is below the range of cosmological relevance. 

We have also shown that previous claims (realized in the two-flavour approximation) 
in the sense that the RG running could provide a substantial enhancement of the 
atmospheric mixing, sin^ 29i cannot work in practice, since the mechanism leads to an 
unacceptably large "solar" splitting, Am2i. We have shown that, unfortunately, this 
is also the case in the more general 3-flavour scenario. 

For the general complex case, most of the conclusions are similar, but there are 
interesting new effects. In particular, the previously considered "radiative" SAMSW 
scenario requires — 0' ~ (unless the "solar" 63 angle is set by hand at a small value 
from the beginning). On the other hand, if — 0' is different from zero and the neutrino 
spectrum has a two- fold degeneracy, the RGBs will generically drive its value to ±7r. 
It is also worth-noticing that the CP phase S can be driven to maximal values by RG 
corrections. On the other hand, if cos 5 7^ 0, the viability of the VO scenario when 
mf ~ 1712 requires, besides the previous conditions, a delicate cancellation involving 
the three CP phases, which makes the scenario more unnatural. 

Addendum 

Shortly after the completion of this work, there appeared a paper by P.H. Chankowski, 



W. Krolikowski and S. Pokorski |^ dealing with similar subjects. More precisely. 



they work out the RGBs for mass eigenvalues and mixing angles below the high energy 
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scale A, assuming that the effective theory below A is the SM or MSSM and a real 
MNS mixing matrix, U. Therefore, their work corresponds to the issues studied here in 
sect. 3. As far as we have checked, their results are in agreement with ours. They also 
notice and stress the existence of stable (infrared pseudo-fixed points) in the evolution 
of the mixings. 

There is a point however in which we disagree. They work out the 2-flavour ap- 
proximation, reaching an RGE for the (atmospheric) angle 6i which coincides with our 
equation (|80|) [or equivalently eq. (^1])], but they argue that the equation normally used 



in the literature, eq. (0), is incorrect and inconsistent with (|80| , |81|) , which is not the 
case (once the difference between mass matrix entries and mass eigenvalues is taken 
into account). 

On the other hand they argue, correctly, that the usual claim that a maximal angle 
is RG stable is not correct. Actually, we have shown that this is also the case in a 
3-fiavour scenario, giving the more general conditions which would guarantee stability. 

Finally, let us remark that in our paper we study also the general complex case, 
i.e. when CP phases are present. In addition, in sect. 2 we work out the RGEs for a 
completely general case (with no reference to e.g. SM or MSSM), showing explicitly 
the appearance of stable (infrared pseudo-fixed points) forms for U. 

Appendix 

We derive here the general form of the stable (infrared pseudo-fixed point) MNS matrix, 
U, in the presence of quasi-degenerate neutrino masses. The generic RGE of the 
neutrino mass matrix is given by {t = log/i) 

^:^ = -(nuM, + M,P + P'^M,), (120) 

from which the RGEs for the mass eigenvalues rui and the mixing matrix U are obtained 
as 

—-1 = -2miPii - miReinu), (121) 

and 

^ = UT. (122) 

dt ^ ^ 



T is an anti-hermitian matrix given by 
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Tij = VijRe{Pij) + i[Vij] ^Im{Pij) + iQij, i 7^ j, 



(123) 



with P = \U^{P + P^)U, Q = -if/t(P - Pt)f/ and V 



V 



\m, 



+ mj)/{mi — rrij). 



Let us suppose that two mass eigenvalues, say mi and 1712, are almost degenerate 
at the initial high-energy scale A. Then, near the starting point of the running, t ~ to, 
T is dominated by the (real) terms proportional to V12 



/ 

T ~ ReiT^i] 
\ 

Then, the formal solution to eq. ( |122| ) is given by 
U' = U{t) ~ ?7exp 




(124) 



gRe{T2{)dt 

SlRe{T2i)dt 





t/i?i2(r), (125) 



where i?i2(r) is an ordinary rotation in the 1-2 plane by an angle F 

Ri2{T) = 




(126) 



The value of T will be such that dU/dt becomes non-singular, i.e. it renders Re{Ti2) 
0. Thus the stable matrix U' satisfies 



Re[P[2) = ^Re{U'^{P + P^)U')u = , 
which we can solve for F in terms of initial quantities obtaining 

tan2F '""^^^^ 



(127) 



'128) 



P22 — Pu 

Therefore, using (|125|) , we get unambiguously the stable U' matrix in terms of the 
initial U. 

One can further justify the condition ( |127|) by showing that the RG evolution actu- 
ally drives Re{Pi2) towards zero, as expected. The relevant RG when V12 dominates 
is 

j^R<Pi2) - (Ai - p22)VuRe{Pu)- (129) 

For a sufficiently long running interval one has [using (|121|) ] 

1 



V 



12 



(Pn-P22)log(A/^) 
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(130) 



and inserting this in ( |129| ) we get Re{Pi2) ^ in the infrared. 

Let us consider now the only remaining possibihty, i.e. that all the mass eigenvalues 
mi, 1712, 'mz are (almost) degenerate. Similarly to the previous case, for t r^ to, T is 
dominated by the (real) terms 



T 




So, the solution to eq. 
U' = U{t) ~ [/exp 



-Re{T2i) 


i?e(T32) 

is formally given by 






(131) 



tRe{T2i)dt 
jlRe{T,,)dt 



gRe{T2i)dt 


-JlReiTs2)dt 



Jl ReiTs,)dt' 

gRe{Ts2)dt 





where i? is a general 3x3 rotation, depending on three angles, say Fi, r2, T^, 



C2C3 
R = I -C1S3 - S1S2C3 

■51^3 — C1S2C3 



C2S3 S2 

C1C3 — S1S2S3 S1C2 
— S1C3 — C1S2S3 C1C2 



UR, 

(132) 

(133) 



with Si = sinFj, q = cosFj (the resemblance with the 'CKM' matrix is only formal). 
R is determined, as in the case of twofold degeneracy considered before, by the fact 



that dU/dt becomes non-singular, i.e. Re{Tij 

1 



Re{P[ 



VI) 



0. Thus 

i?e(t/'t(P + P^)t/')i2 = 0, 



1 



i?e(P^i) = -ReiU'^iP + P^)U%, = 0, 
^e(P^2) = lReiU'\P + Pt)f/')32 = 0. 



(134) 

(135) 
(136) 



Substituting U' = UR, these equations determine completely the rotation R, and 
therefore U', in terms of U and the initial parameters. Another interesting remark is 
that if one has information on the mixing angles in the stable matrix U', the phases 
will be constrained, so as to satisfy the set of equations ( p.34| - |136|) . And viceversa, if 
the phases are known, the mixing angles cannot take arbitrary values. 

Again, one can show that in this case RG evolution drives Pe(Pjj) —>■ for all i,j. 
The relevant RGE that shows this is now 

~ {[Pe(A2)]' + [Re{P23)? + [Re{P3i)]'} = (Ai - P22)Vi2[Pe(A2)]' 

+ (A2 - p33)^23[Re{p23)? + (As " Ai) V31 [Pe( Ai)]'- (137) 



2dt 
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Using 

Vii ~ —. . , (138) 

'' (P,,-P,,)log(A//x)' 

in the equation above, we get Re{Pij) ^ in the infrared for all i,j. 
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